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I. INTRODUCTION 



The principle of thermodynamic similarity is one of the fundamental guiding principle in 
investigation the thermodynamic properties of matter. It allows to determine the classes of 
the substances with similar thermodynamic behavior based on the similarity of the inter- 
particle interactions, which in its turn is based on the similarity of the particles (electronic 
shells, molecular excitation spectra etc.) and as a sequence the similarity of the interparticle 
interactions. The existence of the universality classes of critical behavior can be considered 
as the extension of the thermodynamic similarity principle to the long range fluctuations 
since the short range effects does not change the asymptotic behavior of the thermodynamic 
quantities. In particular, for the systems obeying the the thermodynamic similarity principle 
the dimensionless critical amphtudes should coincide. The demand of the thermodynamic 
similarity and as the consequence the principle of corresponding states (PCS) is very strong 
condition which in fact assumes the similarity of the energy spectrum of the systems since 
the thermodynamics depends on both the external and the internal molecular degrees of 
freedom. 

It is well known that the thermodynamics is determined by the integral characteristics 
of the interaction potential. The classical PCS states the trivial fact that the systems with 
if the interaction potential is of two-particle character and has the form: 



where Uq is the depth of the potential well and tq is the characteristic molecular size corre- 
spondingly. Then the corresponding equation of state (EOS) has the universal dimensionless 
form 



where p* — Ptq/Uq is the dimensionless pressure P, T* — T/Uq and v* — v/vq are the 
dimensionless temperature T and the specific volume v correspondingly. The famous van 
der Waals (vdW) equation of state is the illustration of the PCS. Throughout the paper we 
use the units where ks = 1- In particular the locus of the critical point (CP) depends on 
the potential depth Uq and the molecular radius (hard core radius) and it is expected that 
the locus of the triple points are the same too. 

Sure that such a simple form of the interparticle interaction is nothing than oversim- 
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TABLE I: Basic geometric and energetic parameters of the dimers and the molecules. Notations: 



^vdW ; 



is the van der Waals radius of the molecule, hond bond lentght of the dimer, Sdiss/Tc is the 
dissociation energy of the dimer in units of Tc, is the compressibility factor. 

plification since the classical potential Vir) is the effective interaction. It is obtained via 
integrating out the electronic degrees of freedom and in fact becomes dependent on the 
atomic number A and the de Boer parameter A = ^.^J^^^ (note that m = Anip, where nip 
is the proton mass). For the noble fluids the latter varies in wide range from A^e = 2.7 to 
Axe = 0.06 i i i. 

Nevertheless the quantum effects manifest themselves in dual way: a) the corrections to 
the asymptotic of classical interaction 0]; b) the formation of bound states (e.g. dimers, 
trimers etc.). In the flrst case the potential conserves its structure in essential. The bound 
states lead to the additional contributions to the basic thermodynamic quantities (the spe- 
ciflc volume, heat capacity etc.). The system becomes the mixture of monomers, dimers 
etc. The formation of dimers in noble gases is known [Q. They manifest itself in various 
effects based on the anisotropy structure of the dimers, e.g. in Raman spectra 0, §, 



depolarization of the light scattering ||^, 10 1 etc. But their existence in liquid phase near the 
CP has not been discussed. Note that the dimers have flnite lifetimes but in equilibrium 
there is certain amount of them. 

As it follows from the Table | the number of the vibrational levels Nmb{Uo) is different 
for heavy noble gases. The dimerization breaks the PCS. The manifestation of this effect is 
expected in gaseous phase of low density (see below). But the degree of dimerization is low 
for dense liquid states. Only in the vicinity of the critical point the contribution from the 
dimers could be great enough to break the apphcability of the PCS. 



3 



Therefore the equation of state (EOS) takes the form: 

p* = p{T*, vl, v;,...)+ Pg{T*,A, vl, v;,...), (2) 

where subscript "q" stands for quantum correction and v* are the dimensionless specific 
volumes of monomers, dimers, etc. 

Because of the complete electronic shells and as a sequence the spherical symmetry of 
the ground state the noble fluids are ideal substances for the test of the PCS. 

The aim of this paper is the analysis of the applicability of the PCS for the noble fluids 
and the interpretation of the data from the point view of dimerization processes in these 
substances. More exactly we analyze 1) the applicability of the PCS for the description 
of the noble gases on their coexistence curves; 2) the manifestation of the dimerization 
effects on the temperature behavior of their the rectilinear diameters; 3) the influence of the 
dimerization on the peculiarities of the critical fluctuations. 

In accordance with this Section || is devoted to the careful analysis of the experimental 
data (all the data are from NIST open database). In Section |IT| the purposeful choice of 
the reference system is discussed. In Section |V] we consider simple model to account the 
dimerization process for noble fluids. The mean field behavior of the rectilinear diameters 
are considered in Section |IV| and the critical anomaly is discussed in Section |V1[ 

II. EXPERIMENTAL FACTS 

In this Section we discuss: 1) the comparable behavior of the coexistence curves (CC) 
for noble fluids ; 2) the behavior of the number density along the CC. Further we will use 
the variables reduced to the CP values: 

p = P/Pc, f = T/Tc, h = n/nc. 

First we consider the role of the quantum effects in the EOS. On Fig. we give the CC in 
"p— T" coordinates for classical the noble fluids (A^e, Ar, Kr, Xe) and the quantum ones {He 
and H2). We see that in "p — T" coordinates the classical noble fluids obey the PCS. The 
deviation from the classical behavior is connected with the quantum effects in accordance 
with Eq. d^). They are greater for He than for molecular hydrogen. Note that "p — T" 
diagram is the robust thermodynamic characteristics which does not allow to distinguish 
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FIG. 1: The CC in p — T coordinates. 

fine details of tlie tliermodynamic beliavior. Tlie variables fi, T are more suitable in tliis 
sense. 

Let us consider tlie dimensionless ratio R{2,1) = n2/ni of the dimensionless densities 
of noble fluid "2" to the density of noble fluid "1", which in some reason is chosen as the 
reference fluid. The results of comparison for the liquid and gaseous branches of Kr and 
Xe with respect to Ar are presented in Fig. 0. As is seen from it the deviation from unity 
for liquid branch within the range of the experimental errors (< 1%). For the gaseous phase 
the difference between the Kr and Xe curves reaches 3%. 

The deviation from unity can be considered as the measure of the violation of the PCS. 
In accordance with said above the deviation is noticeable in the gaseous phase where the 
dimerization is more explicit. The physical ground for the existence of the dimers in near 
critical region is the almost free rotation of the dimers allowed at low enough density. As is 
known this leads to the existence of the metastable states for the dimers with a wide range 
of life time because of the centrifugal barrier |l2l. In liquid phase the deviation becomes 
essential in the near critical region. In order to clarify the physical nature of such deviation 
we take into account that along the CC the density is the sum of classical and quantum 
contributions. According to Eq. (0) for the density we have: 

n = nd{f)+fi,{f-K) (3) 

Note that the quantum correction nq{T; A) is monotonous decreasing function on T. As is 
follows from the results presented in Fig. |I] for classical noble fluids the quantum corrections 
are negligibly small. Therefore the deviation R from 1 in liquid and gaseous phase is naturally 
to connect with the dimerization. In order to estimate the dimerization degree it is useful 
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FIG. 2: The ratio R{Kr, Ar) (dashed) and R{Xe, Ar) (bold) 



to search for the reference fluid in which the dimerization is depressed. 

In the following section we show that molecular hydrogen or A^e can be chosen that such 
fluid. 



III. REFERENCE SYSTEM 

Let us give the arguments that 1) the molecular hydrogen H2 is the simplest semiquantum 
fluid and 2) A^e is the simplest classical fluid. This means that in liquid and gaseous phases of 
these fluids in near critical region the amount of the dimers is very small and does not affect 
the thermodynamic behavior. It becomes clear from the value of the dissociation energy 
Ediss reduced to the critical temperature Tg. In molecular hydrogen the dissociation energy 
for its dimer is an order of magnitude smaller than the critical temperature eaiss/Tc ~ 0.11. 
In He no dimers exists except for the giant states [|T3| which is out of the scope of the paper. 
For Ne the dimerization is weak in comparison to other classical noble fluids. A^e is the 
closest to H2 among the classical noble fluids. The comparison assumes that the classical 
constituent of the density for A^e and H2 is different from that of heavy noble gases. Note 
that their critical temperatures is pretty close ~ AA.5 K and ~ 33. K correspondingly 
and high enough so that the difference of these light noble gases from the heavy ones can 
not be attributed to quantum corrections which negligible in near critical region. 

In accordance to Table | the formation of the dimers takes place mainly in classical noble 
fluids. Though the formation of dimers for A^e is possible the number of vibrational levels 
is not big and the dimerization can be neglected in comparison with other classical noble 
gases. 
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FIG. 3: Residual specific heat (in units oi ks, Cy = 3/2 is the specific heat of the monoatomic 
ideal gas) of the noble gases along liquid (black) and vapor (gray) branches of the binodal. The 
inset shows the increment of the specific heat caused by the dimerization in near critical liquid 
phase (bold segment of the curve). 



In order to estimate the fraction of dimers in liquid phase we use the specific (per atom) 
heat data (see Fig. §). Here we see the monotonic change of the residual (after subtraction of 
the specific heat of the monoatomic ideal gas contribution) specific heat in gaseous phase due 
to increase of interparticle interaction with the density up to the temperature T/Tc ~ 0.95, 
which has the meaning of the Ginzburg temperature Jl^. Note that in this phase for 
T/Tc < 0.95 the specific heat of Ar is noticeably less than that for Kr and Xe due to more 
strong interpaticle interaction. The opposite behavior of the specific heat takes place in liquid 
phase. Since the density of the liquid phase decreases along the binodal the corresponding 
decrease of the residual specific heat is observed in the temperature region T/Tc < 0.95. For 
T/Tc > 0.95 the specific heat of both phases diverges due to critical fluctuations. 

From the physical point the decrease of the specific heat of the liquid phase outside of the 
fluctuational region is connected with the decrease of the average interparticle interaction 
as the density decreases. The existence of the minimum at T^m ~ 0.87 and the subsequent 
mcrease of Cv until Tdiss ~ 0.95 (see Fig. D can be attributed to the formation of the 
dimers. Representing the specific heat as 



Q 7 

-Ni + -N2 
2 2 



(4) 



where Ni is the number of atoms and N2 is the number of the dimers, we get the estimation 
for the degree of dimerization (the ratio of the number of dimers to the total number of 
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FIG. 4: The difference of the entropies of the hquid phases of Ne and Xe. 
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where A^o = + '^^2 is the total number of atoms, ACy = CyiTdiss) — Cy(Tmin) is 
the increment of the specific heat due to the formation of the dimers. Here it is assumed 
that the dimer has 6 degrees of freedom (3 translational, 2 rotational and one vibrational). 
Of course the dimerization in liquid phase should manifest itself in the behavior of the 
entropy too (see Fig. ^). Here we see that up to the temperature T 0.87 the monotonous 
behavior of Snc — Sxe takes place. It signifies about the absence of the clusterization 
ordering in Xe. Noticeable diminishing the monotonic decrease in the temperature interval 
0.87 < T/Tc < 0.95 is naturally explained by the dimerization of Xe leading to decrease of 
its entropy in comparison with that of A^e. 

This energetic condition should complemented with the geometric condition of free rota- 
tion of the dimers. Two molecules occupy in the average the volume 2/n, therefore he radius 

1 /3 

of the effective cell for a dimer is r^ff — {^^^ ■ The rotating dimer occupies the volume 



of the radius = The ratio Teff/rd which is shown in Fig. y and becomes greater than 
unity in the vicinity of the CP. 

Let us consider one of the light noble fluid as the reference system. One can conclude 
that the ratio of the densities for the heavy and light noble gases: 



is monotonous along the coexistence curve provided that the classical contributions are the 
same. In Fig. ^ we give the ratio R for heavy noble fluids {Ar, Kr, Xe) with respect to the 



(5) 
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lighter ones {H2, Ne, He). It shows that in the near critical region according to the molecular 
hydrogen is more like A^e than He. The additional support in favor of the thermodynamic 
similarity H2 and A^e is the closeness of their CP locus. Note that in gaseous phase the 
significant deviation of R from 1 with lowering the temperature (density of saturated vapor 
decreases too) is connected with the dimerization in accordance with what has been said 
in the Introduction. From Fig. |TT| we can conclude that the dimerization degree in gaseous 
phase is more than 10%. 

From what has been said above the nonmonotonic change of R can be attributed to the 
temperature behavior of the classical part Vd of the specific volume only. The classical be- 
havior for the noble fiuids is not the same as for H2 because of formation of the dimers. That 
mean that the classical part of the thermodynamic functions corresponds to the monoatomic 
vdW fiuid. 
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FIG. 5: The ratio i?j along the binodal for Ar (dotted), Kr (dashed), Xe (bold) and different 
reference fluids (see the captions). 



The dissociation energy for dimers of heavy noble fiuids is close to the dimer bound energy 
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FIG. 6: The dimensionless distance ^ = 2 (||- ) / {Ibond + f'") 

(the critical temperature shghtly (4% — 6%) higher than the dissociation energy of the dimer 
^diss/Tc ~ 0.95, see Table |) ||^. From Fig ^ one can see that nonmonotonic change in the 
ratio R occurs right at the temperature kT ~ Sdiss- 

It can be seen that the geometric (see Fig. P) and energetic conditions give the ground 
to expect the nonzeroth dimerization of heavy noble fluids in the vicinity of the CP. One 
can choose either A^e or molecular hydrogen as the reference fluid since they can be treated 
without taking into account of the dimers Ne2 or {H2)2- In although the energetic 
condition allows the formation of dimers the number of the bound states is low and its 
dimerization in near critical region is negligibly small in comparison with heavier noble 
fluids. The spectrum of dimers of heavy noble fluids have many vibrational states. Thus the 
dimers of heavy noble fluid are more stable. From this point of view the molecular hydrogen 
and A^e are indeed the "simplest" molecular liquids. The quantum effect for H2 should be 
taken into account via the representation Eq. (|^). 

From here we can conclude that the formation of dimers in liquid phase takes place in 
near critical region where the density is low enough so that both energetic and geometric 
conditions are valid. Therefore there is the difference in short range correlations in light and 
heavy noble fluids which is based on the formation of bound clusters. 

In fact, the liquid phase of a noble gas is the mixture of monomers and small quantities of 
dimers. Since their concentration depends on the kind of matter, the critical amplitudes for 
different noble gases will take values deviating from those, which follow from the principle 
of corresponding states. In this sense one can say that the quantum properties, namely the 
stability properties of bound states, influence the critical parameters of the classical noble 
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fluids. This result was noticed in numerical simulations 
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IV. THE BEHAVIOR OF THE RECTILINEAR DIAMETER 

The rectilinear diameter for the density: 

is important characteristic of the phase coexistence. Such a quantity measures the asym- 
metry of the phases with respect to the order parameter chosen. This asymmetry is related 
to the particle-hole asymmetry [|16] and is sensitive to the interparticle correlations which 
governs the expansion rate of the phases below the CP. 

Let us make some general comments about the temperature behavior of the rectilinear 
diameter. Obviously one can distinguish between three characteristic regions: a) fluctuation 
region Tqi < T < T^. where Ta is the Ginzburg temperature |T^; b) normal states near the 
triple point T > Ttr] c) the intermediate one Tfr < T < Tc. The temperature dependence 
of the rectilinear diameter of the density is governed by the different physical mechanisms 
in theses regions. Near the triple point the hard core interactions in liquid phase play the 
dominant role so one can expect that the behavior in this region should be determined by the 
hard core interactions only. The near critical behavior is governed by long range fluctuations. 
The specific temperature behavior in the intermediate region depends on the balance of 
the short range and long range correlations. Thus the rectilinear diameter characterizes 
quantitatively to what extent the rate of short range ordering in liquid phase overcomes the 
rate of disorder in coexisting vapor. 

These general remarks are illustrated by the data for the rectilinear diameter of the 
density (see Fig. |^). The quantity pd is sensitive to the order in the local configuration. 
In mono atomic liquids like noble ones the monotonic change of this quantity takes place. 
This is because of the stronger correlations in liquid phase. The role of the long range 
correlations in the formation of non monotonic behavior of Ud in the near critical region is 



discussed below in Section VI 



In the fluids with polyatomic molecules the volume available for the thermal motion of 
the molecule modifies the rotational motion in the effective "cavity". It leads to the fact 
that the contribution of the inner degrees of freedom, in particular the rotational ones, 
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FIG. 7: The rectilinear diameter for the density na- 

becomes dependent on the specific volume. As a consequence the critical amplitudes should 
be different and this is demonstrated by the behavior of the rectilinear diameter in near 
critical region (see Fig. ^). Note that the temperature behavior of the rectilinear diameter 
of the density is the fine criterion for the validity of the EOS. 

The natural characteristic of the ordering process in the coexistent phases is the en- 
tropy. Therefore the qualitative difference in the ordering in coexisting phases can be also 
characterized by the rectilinear diameter for the specific (per particle) entropy: 



Si + Sg 



Sr. 



(7) 



The behavior of S^, as well as ria, for different substances which obey the PCS should be the 
same. For comparison on Figs. ^ we show the data Sd for noble fluids. Except the critical 
region the behavior of Sd within the mean fleld approximation differs from the experimental 
curve by the shift of order AS. The appearance of this difference is naturally explained by 
the the fluctuation shift of the critical temperature: 



AS = -Sr^ ct^^ In 



(™/) 



where ch^^^ is the dimensionless value of the mean fleld value of the speciflc heat, which for 
the heavy noble fluids takes the value c""''^'^ ^2). So that the entropy diameter value at the 
CP is less than that in the mean fleld approximation. 

Let us state the relation between the behavior of diameters of the entropy Sd and the 
density n^. We start from the basic thermodynamical representation of the entropy in 
dimensionless units (in units of the Boltzmann constant ks)'. 



T 

S = Sc + c^ln— + f{n) - f{nc) 

J- r. 



(8) 
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(a) Ar, Kr, Xe (b) He, H2,Ne 

FIG. 8: The rectilinear diameter for the entropy S^- 



where Sc is the entropy at the CP, c„ is the dimensionless specific heat, f{n) is some function 
describing the density dependence of the entropy. From Eq. (R|) one can conclude that: 



Sd 



SI) , Sa) 



2 T, 2 

In the mean filed region below the CP we have: 



In — + 



(9) 



T — T 

±6A/|r[ + a|r| + . . . , r = < 

J- r 



(10) 



where 6n = n — ric, and the indices l,g stand for the coexisting liquid and gas phases 
correspondingly. The searched relation between Sd and Ud near the CP is obtained by the 
expansion of Eq. (y) with respect to the density: 



Sd = n, 



dS 
On 



n-d + 



m — + 



T ' 4 \ dn^ 



6nf + 



dm? 



6ni 
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Further we use the thermodynamic relations [|l^ 



d_S_ 
dn 
d^S 



1 dp 
2 dp 



< 0, 



d'^p 



rii dTdn 



(12) 
(13) 



Omitting the irrelevant terms near the critical point we get: 



Sd = - 



1 dp 
n'c df 



{a-b')+b' 



d'^p 
dndT 



+ 



„(0 , J9) 



r 



+ o{\t\). 



(14) 



The coefficient h depends on the parameters of the EOS. In particular for the vdW EOS 
,6 = 2 and ^ = 0.6, for the Bertlo EOS a = |,6 = 2^2 with % = 1.35 in near 
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(a)The behavior of the residual entropy (b)The behavior of the diameters for the 

gases. total and residual entropy 

FIG. 9: The rectilinear diameter Sd of the total entropy and the residual one S^^^^ of the noble 
gases A^e (doted) and Ar,Kr,Xe (bold). 

critical mean-field region. The slope of Sd increases with the density asjTumetry which is 
described by the coefficient a. For lighter noble fluids the slope of the entropy diameter is 
opposite to that for the heavy noble fluids. This is because of the fact that for light noble 
fluids the density asymmetry is weak (see Fig. |^ and the specific heat for the liquid phase 
is smaller in comparison with that for the heavy noble ones (see Fig. |^. 



The difference between the behaviors of 5*^ for light and heave noble gases is explained by 
the dimerization effects (see Fig. |]). The calculation of the diameter for the total entropy 
Sd as well as for its the residual part: 

■,{int) _ + 



where 



niint) _ ~r Q(mi) 



S^''^'^ = S- S^'"'^ , = ^ In ^ - In - 

2 Tr Tlr 



with the help of the vdW EOS, shows that the slope of the 5*^ increases with account of the 
dimers (see remarks about Eq. (|14|)). The comparison of the experimental data for Ne and 
Ar, Kr, Xe supports this conclusion (see Fig. |^) In order to clarify some qualitative details 
of the temperature behavior of 5*^ we will discuss the dimerization of the heavy noble fluids 
in more details in the next Section. The behavior of Sd in the fluctuation region we will 
discuss after this. 
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V. THE DEGREE OF THE DIMERIZATION FROM THE CHEMICAL EQUILIB- 
RIUM 



First let us make some general conclusions about the temperature behavior of Eq. 
along the binodal. For the description of the dimerization let us introduce the quantity 

^dim = — , n = -TT , ni = —- (15) 

rii n('') n^^) 

where n is the total density (dimers and monomers) and rii is the number density for 
monomeric fluid. 

For high enough density in liquid branch for the substances of interest yidim will be close 
to unity far away from the critical point. It is seems to be natural that with increasing the 
density the dimerization becomes negligible. 

At T = Tc according to the definition = 1- The monotonous change of ^idimiT) 

between these boundary values shows the quantitative similarity between substances and 
thus the validity of the PCS. In contrast to this the nonmonotonic behavior for '^y^dimiT) 
along the binodal signifies the qualitative difference in the process of thermal expansion. 
According to what has been said above the dimerization takes place at low enough densities 
near the critical point. The formation of the dimers lead to the flattenning of the binodal 
thus increasing the difference in the densities of the coexisting phases. The increase of 
the density in liquid phase far from the critical point makes the dimerization impossible 
both from the energetic and configurational point of view. In such a case one expects the 
nonmonotonous behavior of yidim(T) along the binodal in liquid phase. 

Here we calculate the dimerization degree for the heavy noble fluids using the chemical 
equilibrium for the dimerization process: 

monomer + monomer ^ dimer 

Let us introduce the degree of association a: 

ni = {l-a)no, ris = | . (16) 

where n2 is the number density of the dimers and no = ni + 2n2 is the total number density 
of atoms. The equation of chemical equilibrium for a is: 

(n (ex) (ex) , 7-1 \ 
f ) 
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where Ediss the dissociation energy of the dimer (in calculation Ediss = 0.97 Tc. The excessive 
chemical potentials are: 



(ex) 



drii 



where F'^^^^ is the excessive free energy is the difference of the total free energy F and the 
ideal part F'^'"'^: 

p{ex) ^ p _ p{id) _ ^-j^g^ 

In its turn F^^^^ can be decomposed 

p{ex) _ p{hc) _|_ piint) ^ 

where F^^^^ is the hard core contribution and F^**^*^ is the contribution due to van der Waals 
interaction. Further we will use the quantities reduced to the vdW values of the critical 
parameters for monomeric fluid. The interaction term F*^*"*-' we take in standard binary 
form: 
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piint) ^ _ (^2^^^^^2+7^2) 



(19) 



Assuming the similarity of the interparticle interactions (see Eq. (P) for monomers and 
dimers, within the vdW approximation and A and 7 are determined as following: 



A 



U 



12 



1 + ^ 



7 



f^22 / Cr2 



8f/ii V 0-1/ Uu \ai 

where f/n, ?7i25 f^22 are the interaction constants for " monomer- monomer" , "monomer- 
dimer" and " dimer- dimer" correspondingly. In the order of magnitude A > 1 and 7 > 2. 
The hard core contribution is taken in vdW form with account of the mutual influence of 
excluded volume for the monomers and the dimers: 



p{hc) 



-T 



rii 61 \ 1 / , b2n2 hi rii 
n, ln( 1---^ +n2 In 1 - ^ ^ 



(20) 



Here hi and 62 are the specific volumes per dimer (in units of the monomer volume) in 
environment of the monomers and the dimers correspondingly. As the calculations show that 
the account of the dimerization leads to the broadening of the binodal in n — T coordinates 



(see Fig. 10) 



The ratio 91 is shown in Fig. |lT|. The comparison with the experimental data of Rne 
for the pair of substances Xe and A^e as the vdW fluids with dimers and without them 
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FIG. 10: The binodals for vdW EOS (dashed) and for vdW with account of the dimerization (bold) 

correspondingly. It shows that the mean field approach is adequate at least qualitatively 
(see Fig. |Tl|). 

The nonmonotonic behavior of i?7Ve as it follows from Fig. ^ takes place in close vicinity 
of the critical point. The formation of the dimers becomes impossible if the specific volume 
becomes less that 2vo (in reduced coordinates n ^ 1.5). Simple estimate shows that it occurs 
exactly at the temperature T ^ 0.95 ^ 0.97 T^. In the proposed mean field description dldim. 
changes in broader region, h ^ 2 (see Fig. |TD|). Such a discrepancy can be attributed to 
the inadequate account of the hard core effects which mainly governs the destruction of 
the dimers when the density increases. The use of more appropriate approximation for the 
hard sphere contribution, e.g. Carnahan-Starling, decreases the value of yidim beyond the 
fluctuation region. The behavior of in the mean field region, T/T^ < 0.9, depends on 
the short range correlations. Obviously the difference in the character of such short range 
correlations is responsible for the different behavior of the rectilinear diameter of the entropy 
Sd which is clearly shown on Fig. ^ for the heavy {Ar, Kr, Xe) and light noble fluids. As 
has been noted above the existence of the dimers is not possible in this region. But the 
stronger short range correlations, as compared to those in lighter ones, should be expected 
in the heavy noble gases as the remnants of the near critical point dimerization. 
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FIG. 11: The comparison of the theoretical (dashed) result for 9\dim and the experimental 
R{Xe, Ne) (bold) data (the parameters 61 = 1, 62 = 3, 7 = 2, A = 1). 

VI. THE RECTILINEAR DIAMETER ANOMALY 

The asymptotic behavior of the rectihnear diameter of the density in the fluctuational 
region is as following: 

ni + Hg 



rid 



2nr 



1 = ItI^/^ + |r|^^" + Di\t\ + o(r) , Di < . 



(21) 



It contains |rp^-term |T9| , |20| , |2T| . The ground for introduction of such a term as the sequence 
of the choice of the order parameter predicted firstly within the canonical formalism approach 



2^ . Note that according to the ratio of the leading singularity amplitudes for the 



temperature dependence of the rectilinear diameter is "almost" universal: 

Di^ ^ ^1^.-2, < . 

where a is the nonuniversal factor. The quantities fg and Qs are universal functions: 

df six) 



(22) 



fs{x) = fs{-x) , gs{x) 



dx 



They are determined by the asymptotic thermodynamic potential 



|2-a 



fs 



A 



/3+7 



+ 



in isomorphic variables Ai - conjugated field-like variable and A2 - temperature-like variable. 

The temperature behavior of the rectilinear diameter strongly depends on the interplay 
between the values of D^. Basing on the experimental data shown in Fig. |^ it seems natural 
to assume that Di > 0. From the physical interpretation of the it is clear that the 
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repulsive terms in pressure decrease tlie absolute value of Di. Since the 2/3 singularity is the 
strongest one the non monotonous behavior of the density rectilinear diameter takes place 
if D213 < 0. In order this behavior manifests clearly there should be -Di_q, < I-D2/3I which 
according to Eq. (521) means that a < 1. The parameter a depends on the choice of the 



initial order parameter. This explains the fact that 2f3 singularity amplitude varies greatly 
depending on the order parameter chosen. This usually is interpreted as the division of the 



order parameters into more and less symmetrical ones, e.g. like density p and volume v p4 

Even if > the 2(3 singularity can be masked by the interplay of linear and "1 — a " 
terms Di ^ |-Di-a| if D2f3/\Di^a\ ^1 (a < 1, see Eq. (^)). In particular from the results 



of [25] it follows that — —0.1 with > 0. Note that in |25] the data for the liquids 
with different degree of the binodal were analyzed, namely SFq as typical molecular liquid 
and alkali metals Cs, Rb. 

If D213 > and in addition (1 — > Di then essentially nonmonotonous behavior 

for the rectilinear diameter takes place. There is minimum at 

1 , (l-a)Z^i_, 

T • r^i — n 

a Di 

and the maximum 



2(3 D2 



1-2/3- 



1 — a 

The location of this maximum is very close to Tc ( \D2i3/Di^a\ < 1) and can be measured 
only if (1 — a)|Di_Q,| ~ Di. The clear example of such situation is given by rectilinear 
diameter for the binodal of Hg . The data used for light nobel fluids here does not allow 
to distinguish between this case and the situation with < 0. 

The values of the amplitudes D2/3 and -Di-a are the same for the substances which belong 
to the same class of corresponding states (eg. Kr and Xe). As is shown in D2/3 and 
-Di-Q have opposite signs. The coefficient D2/3 can be either positive or negative depending 
on the details of the intermolecular interactions. Taking into account the results of the 



previous section we can conclude that for molecular liquids as well as in liquid metals ||26 
where the clusterization takes place in near critical region the coefficient -D2/3 is positive. 
The same is true for the entropy diameter and any other order parameter for which the 
Landau-Ginsburg functional contains the odd power terms. 

The situation with the rectilinear diameter for lighter noble fluids suggests two possibil- 
ities. The first one is that the amplitude |-Di_a| > (^)-D2/3 > but (1 — a)\Di_a\ > -Di 
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and the binodal is quite symmetrical. The second is that a < 1 so that Di_q,(^) < — -D2/3. 
The data for A^e also conform with such a proposition. Note that helium has "almost" 
symmetrical binodal (|-Di| ^ 1) in comparison with other noble liquids. Obviously this can 
be attributed to the quantum corrections to the equation of state which gives the additional 
repulsive contribution to the pressure. Thus one can expect that in He the coefficient -D2/3- 
From here we can make the following conclusion concerning the sign of the amplitude -D2/3- 
The dimerization influences the character of the short range correlations. If the dimeriza- 
tion is negligible then one can expect that 1^2/3 < 0. For the systems where the dimerization 
takes place D2/3 > 0. This leads to the qualitative difference in the behavior of the rectilinear 
diameter observed. 

Here we perform the analysis of the critical asymptotic of the rectilinear diameter for the 
entropy 

Sd = ^2/3 Irp'^ + |r|i-" + S^\t\ + o(r) . (23) 

The first contribution is caused by the corresponding term of the density diameter and the 
quadratic density contribution in Eq. (|l^). Note the term oc r^~° is generated by the density 
diameter and the specific heat terms according to Eq. (pT]). Note, that as follows from the 
representation (Q) Si^a < 0. In contrast to the coefficients Di^a which is generated by the 
asymmetry of the Hamiltonian, the coefficient Si^a does not vanish in symmetrical case, 
e.g. in the Ising model, and correspondingly: 

= ^JlT^ + Str^ . (24) 

It is obvious that the symmetrical part S'j!!'^'' is generated by the heat capacity singularity. 

As follows from Eq. (|TT]), ^i-q, < 0. In order to distinguish between the r^^ and r^~" 
terms for rectilinear diameter we subtracted the analytic background terms outside the 
critical region. The coefficients Sk of the polynomial 

Sr = Y.s,{TjT-lf 

k=l 

have been determined with the help of the experimental data in the region T/Tc < 0.75. 
The results of calculation of the ratio 82(3 / Si^a are placed in Table H and they confirm 
the conclusion of the canonical formalism about its universal character. Small variations 
in this ratio in dependence on the specific fiuid are explained by the difference between 
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H2 


Ne 


Ar 


Kr 


Xe 


S2l3 
Sl-a 


-0.28 


-0.30 


-0.31 


-0.34 


-0.33 



TABLE 11: The ratio 4^ for the noble gases. 

canonical temperature variable A2 and the reduced temperature r. Besides we neglected 
the symmetrical in S'fj^^^ since it is proportional to the fluctuation correction of the heat 
capacity. The results allow to conclude that S'^j^lf' is smaller that the term generated by 
the asymmetry of the Hamiltonian. 



VII. DISCUSSION 



In this work it is stressed that the leading role in the understanding of the difference and 
similarity in the thermodynamic behavior of fluids belongs to the PCS. The main result of 
the present work is the clear demonstration of that fact, that noble gases do not rigorously 
follow the PCS principle of corresponding states 

As was shown from the analysis of the density and the entropy the deviations from PCS 
are observed 1) for liquid phase near the critical point and 2) and for gaseous phase in 
the whole region of their existence. The latter is especially important circumstance, which 
directly indicates the main cause of such deviations. It is the dimerization processes in noble 
gases. The most essentially the dimerization manifests itself in gaseous phase far away from 
the critical point, where the doubled fraction volume per molecule 2v is noticeably more than 
the own volume for a dimer. The degree of the dimerization here reaches 10%. In liquid phase 
due to the geometrical restriction for the existence of the dimers the temperature interval 
where the dimerization effects in nobel fluids take place is rather narrow due to the closeness 
of the dissociation temperature Tdiss to the critical one Tc {Tdiss ~ 0.96 Tc). Nevertheless 
they are quite noticeable in the temperature behavior of the entropy and the specific heat. 
The analysis of the experimental data with taking into account of the dimerization effect for 
the liquid phase in near critical region give the estimate < 5% for the degree of dimerization. 
Note that these effects can not be adequately described with the standard two parametric 
vdW-like EOS. 

The density and the entropy diameters are the fine characteristics of the EOS. From 
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our analysis it follows that widely used vdW and Bertlo EOS are not able to describe the 
formation of the plato in the temperature dependence of 5*^ for the heavy noble fluids. At 
the same time the density binodal is reproduced quite correctly. 

The analogous situation is characteristic for liquids with nonspherical molecules, first of 
all, for benzene and water. In these cases the behavior of the entropy diameter is essentially 
nonmonotonic. This circumstance for water is expected since the association processes in it 
are considerably more expressed than in noble gases. In ||2^ it was shown that the degree 
of the dimerization near the critical point is close to unit (a ~ 0.9). At removing from the 
critical point the associates in liquid water becomes more complicated. So their influence on 
the flne details of density and the entropy of water is considerably more expressed. Strong 
correlations in liquid benzene lead also to considerable changes in the behavior of S^- The 
speciflc structure manifest itself also in kinetic properties such as viscosity. 

As shown in the work, the dimerization radically influences the behavior of the density 
rectilinear diameter. From the analysis of the last for He, H2 and Ne it follows that 
is negative or takes the small positive value, essentially less than -Di_a. Note, that the 



using of the experimental data from [25] allows to get the following estimate for their ratio: 
D2/3/D,,^ = -0.1. 

It is necessary to emphasize that in accordance with the canonical formalism reasons the 
signs of D213 and -Di_q, should be opposite. This very important conclusion could not be got 
with the help of theory developed in ETI, This question will be discussed in separate 



paper. 

One of the author (V.L.K.) acknowledge the partial flnancial support of the INTAS 
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